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Abstract 

The cross sections of the '^^Tm(Q','>')''-'Lu and '®Tm(Qf,n)'^^Lu reactions have been mea- 
sured first time using a new method, by detecting the characteristic X-ray radiation following 
the electron capture-decay of '^^ '^^Lu. Despite the relatively long half-life of the reaction prod- 
ucts (Ti/2 - 500 and 6.7 days, respectively) it was possible to measure the cross section of the 
'^^Tm(Q',7) '^^Lu reaction close to the Gamow window (Tg = 3.5 GK), between E(, „, = 13.16 and 
17.08 MeV. The '®Tm(a;,n)'^^Lu reaction cross section was measured from Ec.m. - 11-21 MeV 
up to Ec.m. = 17.08 MeV. The experimental results have been compared to theoretical predictions. 

Keywords: Nuclear astrophysics. Nucleosynthesis, Astrophysical y-process. Characteristic 
X-ray emission. Statistical model 



1. Introduction 

In stellar nucleosynthesis involving charged particles two probability distributions are play- 
ing crucial role: the Coulomb penetration probability and the Maxwell-Boltzmann energy distri- 
bution. The first strongly increases with increasing particle energy, but at a given temperature the 
probability of a particle having high energy falls off rapidly, following the Maxwell-Boltzmann 
distribution. The Coulomb penetration probability folded with the Maxwell-Boltzmann distri- 
bution forms the so-called Gamow peak, and the nuclear reactions are occurring in this narrow 
energy window Il]|2]. In order to provide reliable nuclear data for the astrophysical calculations, 
the experiments have to be performed in the Gamow window or at energies as close as possi- 
ble. At higher burning temperatures this energy window can be shifted up to higher energies. 



'Present Address: Laboratori Nazionali del Sud, INFN, Catania, Italia 
Email address: ggkiss®a.tomk.±.im. Fax: (+36) (52) 416 181 (G. G. Kiss ) 



/Physics Letters B 00 (2010) 7-[9] 



2 



Table 1 : Existing experimental (a,y) database relevant for the astrophysical y-process measured via activation method 
and y-counting in the A > 100 region. The target nuclei, the half-lives of the reaction products, the Gamow windows (21 
and the lowest energies measured (En,i„) are given. 



Target 


Half- 


Gamow window 


F 


Ref. 


nucleus 


life [d] 


[MeV] (T9 = 3.5 GK) 


[MeV] 






0.17 


6.6-10.1 


7.57 


m 


"2Sn 


0.10 


6.8-10.7 


9.97 




"^In 


0.12 


6.5 - 9.4 


8.66 


m 


i"Sn 


154 


6.1-9.6 


11.9 




Au 


3.05 


8.5 - 11.7 


17.5 


m 




500.5 


8.2- 11.3 


13.2 


this work 



however, it remains mostly well below the Coulomb barrier, making it hard or even impossi- 
ble to perform experiments at astrophysically relevant energies OJ. Here - since experimental 
data is missing in the most relevant mass range for the so-called astrophysical y-process, too 
- we propose a new experimental approach to measure the cross sections of important nuclear 
reactions. 

In order to understand the synthesis of the so-called p nuclei — the rare proton rich iso- 
topes between ^^Se and '^^Hg |[3| — experimental data at sub-Coulomb energies are particularly 
important. Despite the tremendous experimental and theoretical efforts of recent years, the syn- 
thesis of the p nuclei is still one of the least known processes of nucleosynthesis. It is generally 
accepted, however, that the main stellar mechanism synthesizing these nuclei — the so called 
y-process — is initiated by ('x,n) photodisintegration reactions on preexisting more neutron-rich 
seed nuclei. As the neutron separation energy increases along the iy,n) path towards more neu- 
tron deficient isotopes, ('x,p) or {y,a) reactions become faster and process the material towards 
lower masses ||4][5]|6l. Recently, consistent studies of the nucleosynthesis of the p nuclei become 
available, employing theoretical reaction rates in large reaction networks |5, 6|, proving that in 
the case of the production of the heavy p nuclei (140 < A < 200) the reaction flow is strongly 
sensitive to the {■Y,a) photodisintegration rates. If those reaction rates are high, more material 
will contribute to the synthesis of lower mass p isotopes. On the other hand, if the reaction rates 
are lower, processing toward lower atomic numbers is weaker, resulting in a relative enrichment 
in the higher mass p isotopes such as '^^Hf, '**"W, and '^Opt |6|. Consequently, to reproduce the 
path of the y-process, experimental data is highly needed in this mass range. 

2. Experimental approaches used in y-process studies 

In principle, photodisintegration cross sections can be determined directly by photon induced 
reaction studies Q. However, in such an experiment the target nucleus is always in its ground 
state, whereas in stellar environments thermally populated excited states also contribute to the 
reaction rate leading to large corrections of the ground state rate which can only be modeled theo- 
retically m. It has been shown that the influence of thermal population is much less pronounced 
in capture reactions and therefore it is advantageous to measure in the direction of capture and 
convert the measured rate to the rate of the inverse reaction by applying the principle of detailed 
balance ll9l fT0ll . Recently, several (a,y) cross sections around A 100 have been measured using 
flie activation method |[IT][12l[13l[ll[T5l[T6l|T7l (see Table[T). Above the A ^ 100 mass region, 
however, there are practically no experimental (a,y) data below the Coulomb barrier except the 
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Table 2: Decay parameters of the '*'Tm((!',>')'^^Lu and '^'Tm(Q-,n)'^'Lu reaction products 1211 . 



Residual 


Decay 


Half- 


Energy 


Relative intensity 


nucleus 


mode 


life [d] 


[keV] 


per decay [%] 


i'"Lu 


£ 100% 


500 ± 3.65 


51.35 (K„,) 


43.8 ± 1.4 








52.39 (K„,) 


76.3 ±2.4 


i^^Lu 


e 100% 


6.7 ± 0.04 


51.35 (K„,) 


31.5 ±0.9 








52.39 (K„,) 


54.9 ± 1.5 








810.06 


16.6 ±0.9 








900.72 


29.8 ± 1.3 








912.08 


15.3 ±0.7 








1093.63 


63.0 ±3.0 



Sm(a,y) Gd reaction 1 18 1. The results of this latter experiment (using a method being lim- 
ited by the fact that the yield of emitted a-particles had to be measured) showed that there are 
large discrepancies between the experimental results and theoretical predictions IITSl . 

Though the activation method proves to be successful to measure a-induced cross sections 
it has numerous limitations which can be sorted into two groups (physical and technical limita- 
tions): 

(1) In order to have reasonable count rates, reactions producing nuclei with long half-life 
(more than few days) were typically excluded. Furthermore, the product nuclei must decay via 
the emission of at least one high intensity (typically 70-97%) y transition. 

(2) A detector with a good signal-to-noise ratio is necessary for going to low energies. A 
high efficiency detector must be used to compensate the low cross sections, long half-lives or 
unfavorable y-branching ratios. The laboratory background can be reduced using lead shielding 
but the beam-induced background on the impurities of the target and/or the backing will still 
be present. Furthermore, if the reaction product decays via y cascades, the true coincidence 
summing effect has to be taken into account. 

In this work the cross sections of the '^^Tm(Q',y)'^-'Lu and "''^Tm(Q',n)''^Lu reactions have 
been measured, by detecting the characteristic K„i_2 X-ray lines following the electron capture- 
decay of produced Lutetium isotopes. It has to be realized that unstable nuclei above the mass A 
a! 150 on the proton-rich side close to the valley of stability typically decay by electron capture. 
Such an electron capture is usually followed by X-ray emission. The characteristics X-ray count- 
ing approach has numerous advantages: the relative intensities are high since the electron capture 
is followed dominantly by the emission of either a Kq,; or a Kq,2 X-ray (the probability of Auger 
electron or X-ray emission are typically an order of magnitude smaller). Furthermore, after 
the decay of one product nucleus only one K„ characteristic X-ray is emitted and consequently 
the target can be put close to the surface of the detector the summing effect is not present. 

To measure these X-ray lines a so-called LEPS (Low Energy Photon Spectrometer) detector - 
consisting of a thin germanium crystal with large surface and a thin entrance window - was used. 
This detector has high efficiency for the Kq, lines and insensitive to the high energy y-s of the 
background. Consequently - compared to a standard shielded HPGe detector - the observable 
minimum counting rates were orders of magnitude smaller (see Fig. 1). This way, despite the 
relatively long half-life of the reaction product (T1/2 = 500 days) it was possible to measure 
the cross section of the '^^Tm(ff,y)'^''Lu reaction close to the astrophysically important energy 
region. The "'^Tm(a',n)'^^Lu (T1/2 = 6.7 days) reaction was studied between Ec.m. = 1 1.21 MeV 
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Figure 1: Off-line y- (upper panel) and characteristic X-ray (lower panel) spectra normalized to the length of the 
countings, taken after irradiating a '^'Tm target with 15.5 MeV a beam. The y-lines used to measure the cross section of 
the '"Tm(Q',n) reaction and the characteristic X-ray radiation used to determine the '*''Tm(Q', y) reaction cross section 
are marked. 



uptoEc.m. = 17.08 MeV. 

The disadvantage of the X-ray counting is that this approach is not isotope selective, e.g. 
in this case it is impossible to distinguish between the decays of the different lutetium isotopes 
produced via the '*''Tm(Q',7')'^-'Lu and the '^'^Tm(Q',n)'^-Lu reactions. However, by taking the 
advantage of the vastly different half-lives it is possible to overcome this problem. Namely, 
the X-ray counting was carried out at least two times: at first after the end of the irradiation to 
determine the '^^Tm(Q',n) cross section and about 24 weeks later when the '^^Lu activity of the 
targets drastically reduced and the measured X-ray yield belongs solely to the decay of the '^^Lu, 
produced via the ^^'^Tm(a,j) reaction (for details, see below). 



3. Experimental technique 

Here we provide a proof of concept by measuring the {a,y) and (a,n) reactions on '^^Tm 
at low energies. The targets have been made by evaporating metallic thulium onto 2 yum thick 
Al foils. The thickness of the thulium targets varied between 106 yug/cm^ and 356 /ig/cm^, at 
higher energy irradiations thinner and at lower energy irradiations thicker targets were used. 
The absolute number of target atoms and the uniformity were determined by the PIXE method 
(Proton Induced X-ray Emission) ll20l using the Nuclear Microbeam Facility of ATOMKI. A 2 
MeV proton beam - the beamspot was 3 yum x 3 jjm - was scanned over a surface of 500 yum 
X 500 fim at several different positions on the target. The precision of the determination of the 
number of target atoms was better than 3% and the thickness was found to be uniform within 3%. 
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Table 3: Experimental S factors of the 



'^''Tm(Q',n)"2Lu and "^''Tm(Q-,y)'"Lu reactions. 



[MeV] [MeV] 



[iQ^^MeVbarn] [lO^^MeVbarn] 



11.5 11.21+0.057 

11.85 11.55 + 0.059 

12.2 11.90 + 0.061 

12.5 12.19 + 0.062 

13.5 13.16 + 0.066 

14.0 13.66 + 0.069 

15.0 14.63 + 0.075 

15.5 15.12 + 0.077 

16.0 15.61 ± 0.079 

16.5 16.10 + 0.081 

17.0 16.59 ± 0.084 

17.5 17.08 + 0.086 



24.0+ 1.7 329 ±33 

18.6+1.4 184+16 

12.0 + 0.9 52.9 + 5.2 

9.48 + 1.01 31.0 ±3.2 

6.04 + 0.49 16.1 + 1.2 

4.51 +0.48 7.87 + 0.61 

3.27 + 0.29 4.11 +0.33 

2.23 + 0.20 2.40 + 0.18 



65.7 + 6.9 
53.5 + 5.8 

46.5 + 5.3 

36.6 + 2.6 



Weighing was also used to determine the number of target atoms. The agreement between the 
results of the two methods is within 4%. The '®Tm targets have been irradiated with a beams 
from the MGC cyclotron of ATOMKI. The energy of the a beam was between E = 11.5 MeV - 
17.5 MeV, covered with 0.5 MeV - 1 MeV steps with beam currents of typically 2 fiA. After the 
beam-defining aperture, the chamber was insulated and secondary electron suppression voltage 
of -300 V was applied at the entrance of the chamber The collected charge in each irradiation 
was between 50 mC and 250 mC. After irradiation, Tu.aitmg = Ih - 3h waiting time was used in 
order to let the short-lived, disturbing activities decay. 

The decay parameters of ''^Lu taken from literature are summarized in Table |2] The 
activity of all samples has been measured with a LEPS detector (GL2015R type, produced by 
Canberra). The characteristic X-ray countings were carried out using two geometries: after the 
irradiation, the targets were put 12 cm away from the detector (far geometry). The countings 
were repeated several times, 1-4 weeks later in close geometry, in these cases the target detector 
distance was 5 cm. The absolute efficiency of the LEPS detector was determined using the 
53. 16 keV line of '^^Ba and the 59.54 keV line of the ^'^^ Am calibrated radioactive sources in far 
geometry where the true coincidence effect is negligible. A linear fit for the measured absolute 
efficiencies were used to determine the absolute efficiencies in the 53 keV - 60 keV energy range. 
Taking into account the elapsed time and the decay between the countings, a conversion factor 
between the two geometries was calculated and used to derive the absolute efficiency in close 
geometry. 

At higher bombarding energies (at and above 15 MeV) the y-activity of the targets has also 
been measured with a 40% relative efficiency HPGe detector, while at 13.5 MeV and 14 MeV a 
100% relative efficiency HPGe detector in Ultra Low Background (ULB) configuration has been 
used. The cross section of the (a,n) reaction obtained from the measurements with the LEPS 
and the HPGe detectors were found to be in agreement within 4%. y-lines belonging to the 
'^^Tm(a, y)'^^Lu reaction were not observable in the y-spectra taken either with the ULB or the 
HPGe detectors. 



^An attempt was also made to measure the low activity of some targets in an ultra low background environment in 
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Figure 2: Experimental Tm(Q;.n) 5 factors (exp) are compared to two theoretical calculations. Both calculations use 
the same code |23| with the same settings except for the optical a+nucleus potential. One calculation was performed 
with the standard potential by L24J (Theory 1 ), indicated by solid line (the inset shows the result divided by factor of 2, see 
text). Alternatively, the recent Friihlich potential 1251 |26I was used (Theory2). Plotted is the original, unrenormalized 
result (dashed line). 



The cross sections of the Tm{a,yy \ji reaction were determined by counting the yield 
of the Kffj and K„, characteristic X-ray lines using the LEPS in close geometry. Figure [T] shows 
typical spectra measured with the HPGe detector (upper panel, t„aitina = 18.47 h) and the LEPS 
(lower panel, twaiting = 168 d) after irradiating the '^''Tm target with a 15.5 MeV a beam. In the 
measured X-ray spectra it is impossible to distinguish between events coming from the decay of 
the different Lu isotopes. However, there is a large difference in their half-lives and therefore 
we can assume that - even if the cross sections were the same - the yield measured after the 
irradiation would be dominated by the decay of '^^Lu. To determine the cross section of the 
'^^Tm(Q',y) reaction, the activity measurement using the LEPS was repeated twice, » 168 days 
and 200 days after irradiating the targets with a beam. During the cooling period of more than 24 
weeks the '''^Lu activity of the targets decreased by a factor of more than 3 x 10^, therefore the 
observed X-ray yield belongs solely to the decay of the '^^'Lu produced only by the '^^Tm(Q',7) 
reaction. The {a,y) cross sections obtained in the 2 measurements are within 3%. The target 
material contained Fe and Ge impurities. The 73 keV and 122 keV background peaks, originating 
from the decay of the reaction products of ^"Ge(Q',n)'^Se and ^"^Fe(Q',n)^^Ni, were used to monitor 
the detector efficiency. No long term variation in the efficiency was found. 



the LNGS underground laboratory in Italy. The analysis of these measurements is still in progress and the results will be 
published elsewhere. 
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Figure 3: Same as Fig.jljfor "'Tm(Q',-)')"^Lu. 



4. Results and discussion 

In astrophysical investigations it is common to quote the astrophysical S factor. The cross 
section, cr(E) and the astrophysical S factor, ^(E) at cm. energy E are related by: 



Eo-{E) 



(1) 



exp{2nrj) ' 

with rj being the Sommerfeld parameter. Several irradiations at the same energies were re- 
peated, the astrophysical S factors shown in Table[3]were derived from the averaged results of the 
irradiations weighted by the statistical uncertainty of the measured values. The effective center 
of mass energies (Ec.m.) in the second column take into account the energy loss of the a beam in 
the target. The quoted uncertainty in the Ecm. values corresponds to the energy stability of the a 
beam and to the uncertainty of the energy loss in the target. The uncertainty of the S factors is 
the quadratic sum of the following partial errors: efficiency of the HPGe detector and LEPS (6% 
and 4%, respectively), number of target atoms (4%), current measurement (3%), uncertainty of 
decay parameters (< 5 %) and counting statistics (0.5 - 7%). 

The experimental S factors have been compared to predictions of the statistical model (Hauser- 
Feshbach) code NON-SMOKER*™ v5.8.1 |23|. This code is a further development of the 
NON-SMOKER statistical model code |22| and includes, among other modifications, updated 
masses and excited states of nuclei. The results underline the importance of measurements at 
energies as close as possible to the astrophysically relevant energy window. The reaction cross 
sections are sensitive to different nuclear properties depending on the projectile energy and reac- 
tion type. It is already known that the standard optical a-i-nucleus potential by |24|, widely used 
in astrophysical calculations, may not be optimal at a energies close to or below the Coulomb 
barrier. This can be seen in Fig. |2] showing S factors for the (a',n) channel. In the measured 
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energy range, the theoretical S factors for this reaction are mainly sensitive to the calculated 
average a widths which, in turn, mainly depend on the optical potential used. The calculation 
using the standard potential (Theory 1) overestimates the cross section by a factor of 2 (this is 
demonstrated on the inset of Fig. |2|i. On the other hand, the calculation using the potential of 
Il25ll26i (Theory 2) is much closer to our experimental values, especially at low energies, without 
any renormalization. This confirms earlier findings that this potential describes low-energy a- 
induced reactions better although it does not reproduce scattering data. This puzzle needs further 
investigation both theoretically and experimentally. 

The '^''Tm(a,y) S factors (Fig.[3|^e not described well with either potential. However, the 
values obtained with the potential of [25 "261 (Theory 2) are again much closer to the experimental 
values than the ones obtained with the standard potential |24J (Theory 1). Contrary to the (a,n) 
reaction, the (a,y) reaction cross sections are not only sensitive to the optical potential but also to 
the calculated y width. Since the (a,n) S factors are described well with the potential of [25 /26), 
we can assume that the remaining deficiencies in the reproduction of the (a,-/) excitation function 
are due to the calculated average y width. This width itself depends on the GDR properties and 
the level density of '^^Lu, which are both experimentally undetermined. However, the (a,-/) 
reaction rate at p process temperatures is insensitive to the y width and therefore an inadequate 
modeling is astrophysically inconsequential. 

5. Summary 

In the present work the cross section of the '*'Tm(Q',y) and '^'Tm(a',n) reactions were de- 
termined well below the Coulomb barrier with high precision (< 10% total uncertainty) in spite 
of the very long half-life (500 and 6.7 d) of the reaction products by counting the characteristic 
X-ray radiation. It has to be emphasized that this method is less sensitive to the usual obsta- 
cles of the activation technique such as low y branching ratios or long half-lives. Moreover, not 
only (a, y) reactions in the region of the heavy p nuclei can be studied via this X-ray radiation 
approach but it can be used also for other applications provided the reaction products decay via 
electron-capture, e.g., in (y,n), (y,p) and {y,a) photoactivation studies. 
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